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Abstract 

We calculate shell-model matrix elements of the axial-charge exchange cur- 
rent operators that have been obtained up to the next-to-leading order from 
heavy-fermion chiral perturbation theory. It is found that loop corrections to 
the soft one-pion-exchange contribution are small (around 10 %) and have no 
significant dependence on the nuclear mass number or on the valence-nucleon 
orbits. These results render further support to the chiral-filtering conjecture. 
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1 Introduction 



The nuclear axial-charge operator, which carries valuable information about the role 
of the chiral symmetry in nuclear medium, has been the subject of a great deal of exper- 
imental and theoretical work (for reviews, see 0,0). The crucial point is that A , the 
time component of the axial current, is predicted to have a large exchange current due 
to a one-pion exchange diagram, the structure and strength of which are dictated by the 
soft-pion theorem Q. According to a general argument by Kubodera, Delorme and Rho 
[||], the impulse-approximation (IA) one-body current j4°(l-body) should receive a large 
meson-exchange-current correction yl°(mec) which is mainly due to the soft-pion-exchange 
two-body current A (soft-7r) : 

A = ^(l-body) + A°(mec) (1) 

where 

A°(l-body) = ^-a-pr ± , (2) 

A°(mec) ~ A (saft-rr) = ( 1 + —) e~ m ~ r (S x + a 2 ) ■ r (r a x f 2 )± (3) 

with r = fx — ?2- Here, p = +P/) with pi and f>f the incoming and outgoing nucleon 
momentum; f n is the pion-decay constant, gA the axial-vector coupling constant, m N the 
nucleon mass and the pion mass. A useful measure of the typical effects of A°(soft-7r) was 
given by Delorme 0]. For a system of one valence nucleon added to the Fermi gas, Delorme 
derived the Hartree-Fock-type effective single-particle operator arising from j4°(soft-7r), and 
evaluated its strength relative to the single-particle operator, eq. (j2j) . This Fermi gas model 
estimate indicates that the original single-particle operator should be enhanced by ~ 54% 
(~ 39%) at nuclear matter density (half nuclear density). Thus, indeed, the two-body 
effect is very large. The replacement of the uncorrelated two-nucleon relative-motion wave 
functions used in the Fermi gas model with typical short-range-correlated wave functions 
would reduce the enhancement to some extent but, even with this reduction, the effect of 
j4°(soft-7r) is still expected to be of a substantial magnitude. 

All the existing experimental information supports the significant enhancement of 
yl°(l-body) due to the soft-pion exchange current (see the extensive bibliographies in |2|). 
The most convincing evidence comes from Warburton et a/.'s systematic analyses || [| of 
the first-forbidden transitions over a wide range of the periodic table. The latest results 
of these analyses § can be summarized in terms of S mec , the ratio of the meson-exchange- 
current contribution to the IA one-body contribution: 

_ {A (mec)) 
<W - ^ (i_ body ))- 
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The best fit to data gives: 

W 16 0) = 0.61(3), 

,5 mec ( 132 Sn) = 0.82(7), 

<W 208 Pb) = 0.79(4). (4) 

These results on $rnec tire certainly compatible with the sizeable contribution of ^4^(soft-7r); 
in fact, the empirically determined enhancements even seem somewhat stronger than an- 
ticipated from the A°(soft-7r) contribution. 

Recent theoretical developments are concerned with: (i) the evaluation of heavy- 
meson exchange diagrams [j?], H|; (ii) the estimation of higher order terms in chiral per- 
turbation H |i~0| ; (iii) possible consequences of the in-medium scaling of particle masses 







All these developments are related to the chiral-nltering conjecture ||, [22| . which 
states that, whenever there is no kinematical suppression, one-soft-pion exchange currents 
should dominate, with other shorter-ranged currents "filtered" off by the nuclear medium. 
Although a reasonably reliable argument for the dominance of A°(soft-7r) in ^4°(2-body) 
was presented in Q, the explicit demonstration of the suppression of other short-ranged 
contributions was left for future studies. The subsequent accumulation of experimental 
evidence in support of the chiral-filtering mechanism has further increased the importance 
of this explicit demonstration. 

Several authors addressed this and related issues by evaluating the Feynman diagrams 
involving heavy- meson exchange [0, |j. The upshot of this approach is that some of the 
short-range exchange currents can be rather large but that they tend to cancel among them- 
selves. The net contribution of the heavy-meson exchange diagrams is to enhance 74°(soft-7r) 
by about the right amount to reproduce <5 mec - Although the idea of heavy-meson-exchange 
currents emerges naturally from the picture of the boson-exchange nuclear interactions and 



from the Dirac phenomenology |T^, 13 1 and although its numerical results are generally en- 
couraging, the convergence property of this approach is a non-trivial question. At present, 
there seems to be no clue to a relevant expansion parameter that controls a multitude 
of complicated diagrams describing short-range exchange currents. Also, the drastic can- 
cellation between the individually large terms tends to make their net contribution quite 
model-dependent, a problem which awaits further investigations. 

A systematic study of 74°(2-body) based on chiral perturbation theory (ChPT) was 
taken up by Rho || with the use of the heavy-fermion formalism (HFF) |l4], |l^, [l6|]. In 
ChPT, the chiral symmetry and well-defined counting rules allow a systematic construction 
of an effective Lagrangian. Furthermore, the use of HFF enables us to carry out a consistent 
perturbation expansion even for systems that involve the nucleons^ 1 . A proof of the soft- 



pion-exchange dominance to leading order of chiral expansion was given by Rho @. The 

#1 From a formal point of view the general application of ChPT to a many-body system involves a subtlety, 
since the system has an extra scale /cf (Fermi momentum) in addition to the usual chiral-symmetry breaking 
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extension of the ChPT calculation to the next order (one-loop order) was carried out by 
Park et al. |1C]. Park et al.'s calculation indicates that the loop correction to the soft-pion- 
exchange contribution is very small, supporting the robustness of the chiral filter mechanism. 
In HFF, baryon momentum-dependent terms appear as higher order interaction terms in 
the chiral counting, and this feature leads to a drastic reduction of the contributions of the 
heavy-meson exchange diagrams. The relation between the heavy-meson-exchange approach 
and ChPT is an interesting problem that warrants further careful studies. 

In pfl , the axial-charge transition matrix elements were calculated with the simple- 
minded Fermi-gas model. To obtain more reliable estimates of the contributions of higher- 
order terms in the chiral expansion, we need to use more realistic nuclear wave functions 
than those of the Fermi gas. Also, according to Towner ||, the matrix elements of the 
heavy-meson exchange currents exhibit significant dependence on the valence orbits. It is 
important to examine whether or not the loop corrections in ChPT shows a similar shell 
dependence. In this note, we calculate the shell-model matrix elements of the axial-charge 
exchange current operators obtained in [1C] up to next-to-the- leading order in ChPT. Our 
results indicate that the loop corrections calculated in the shell model are as small as in- 
dicated by the Fermi gas model estimation. Furthermore, the mass- and state-dependence 
of the loop correction is negligible. These results give further support to the chiral-filtering 
conjecture. In the light of our new results, we shall also discuss the problem of the "ex- 
tra" enhancement in the empirical <5 mec . In particular, we comment on the interrelation 
between the heavy-meson-exchange-current method ||] and the ChPT with in-medium 
mass scaling [O]. 



2 Calculation 

We begin by specifying our conventions and definitions. We introduce, as we did in 
PI , the quantity Xi, which is different from the axial-charge operator A by just an overall 
factor, 

M = ^A° = Mi- hody + M 2 - hody + ---. (5) 
9A 

Using AA is slightly more convenient than using A itself. In eq.((5|), -M n -bod y represents an 

n-body operator, and the ellipsis denotes .M n -body (n > 3), which are ignored. .Mi-body is 

scale A ~1 GeV. However, it seems justifiable to use ChPT in the following limited context. We start with 
the generally accepted paradigm that the diluteness of nuclear matter allows us to concentrate on 1-body and 
2-body responses to external probes, ignoring n-body (n > 3) contributions, and that the effective operators 
describing these responses can be obtained by considering Feynman diagrams involving one nucleon (two 
nucleons) for the 1-body (2-body) operators. Once a one-nucleon or two-nucleon subsystem is isolated from 
the ^4-body system, we use ChPT as a method to calculate the Feynman diagrams pertaining to these 
subsystems systematically and consistently with the basic chiral symmetry of QCD. In the present paper, 
the word "ChPT" should be understood in this context. 
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the impulse approximation term, 



2f 2 

-Mi-body = — o-pr^. 

m N 



(6) 



The two-body exchange current, .M 2 -body> is decomposed into a tree part and a loop cor- 
rection. The tree part is the seagull term considered in the original paper [||. Meanwhile 
according to IfTcf] , the loop correction can be decomposed into a one-pion exchange part (a 
loop correction to the tree part through vertex renormalization) and a two-pion-exchange 
part : 



-A^-body = Mtree + A4i 00 p, 
-Mloop = M ln +M 2 n, 



(7) 



with 



M tree — T ' ^ /tree(^) 

Mi* = T^f^r), 

Ml* = f Ml) (r)+f^ f M2) (r) 

The spin-isospin operators t« and t ( 2 ) are defined as 

Tl X f 2 ) ± f • (<7i + <7 2 ), 



(8) 



r« = ( 

t (2) = (n + r 2 ) ± f • (ai x <j 2 ) 



(9) 



where f is the unit vector of r = r± — r 2 and r is the norm of r. Functions /t ree (^), fin( r 
and / 27r (i,2)W are given in Q, 

1 



/tree(0 

/W(r) 



d_ 

dr 



Airr 



(10) 



c 3 "FT hxee{ r ) 



+ 



f 2 

1 



r/ 



1 + 3^1 



f2-K{%){r) 



167r 2 / 2 dr 

1 d 
167r 2 / 2 dr 

1 d 
167r 2 / 2 dr 



K (r)-^o(r) + (2 + 4^) if 2 (r) - if 2 (r) 



2 5 i^o(r) 



(11) 



(12) 



the constant can be extracted from the isovector Dirac form 



The explicit forms of the functions K{(r) and K{(r) are given in Appendix B of ref . [|T 

As noted in 
factor of the nucleon, i.e., 

f 2 

J 7T 



C 3 f2 



0.04784. 



(13) 
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In shell-model calculations, it is convenient to write the above expressions in terms 
of spherical tensors [17], 



i-Mi-body 
iMi 



V3 — ^ [a, ip 



1(0) 



(14) 



where (m,n) = (2,4) for i = 'tree', l-zr, 2-7r(l), which are proportional to and (m,n) = 
(4,2) for i = 27r(2), which is proportional to T~( 2 \ Here, S„ (t) and Hn\a) are defined in 



(fi x r 2 ) 



4V) 



m + 7- 2 ) 



(15) 



and similarly for En (<r), n = 2,4. The square bracket with superscript (0) in the above 
equations represents the vector coupling of two spherical tensors to a resultant^ 2 ; in this 
case two vectors are coupled to form a tensor of rank zero : 

(o) 1 



o, b 



V3 



a • b. 



(16) 



In the simplest version of the shell model, the ground state (Fock state) of a nucleus 
is taken as a Slater determinant of single-particle states, which here will be chosen as 
eigenstates of a spherical simple harmonic oscillator. In an oscillator Hamiltonian, there 
are two parameters, one is the mass, which will be identified as the nucleon mass, and the 
other is the oscillator frequency, u, which is determined by 



fiuj 



45A- 1/3 - 25A~ 2/3 ) [MeV] 



(17) 



where A is the mass number of the nucleus. For an odd-mass nucleus, the wavefunction is 
written as a simple product of a Fock state for the neighbouring even-mass nucleus and a 
single-particle state, ip a , viz. \if) a ;F). Here the ip a denotes all the quantum numbers of the 
single-particle states indexed by the subscript a, = (n a , l a , s a = jai jaz, ta = \ > taz), 
where n a {= 0, 1, 2, • • •) is the principal quantum number, l a the orbital angular momentum, 
s a the spin, j a the total spin, and the subscript z identifies the third component. The use of 
such a shell- model calculation is given in detail in [17, here we just explain the definition 
of reduced matrix elements and some features of our calculation. 

Consider an operator with definite multipolarity J and isospin T with third compo- 
nent J z and T z , Xj^'^; our axial-charge operator can be regarded as Our goal is 
#2 In general, 

[y(*i)(l), yfe)(2)]^= J2 {kim 1 ,k 2 m 2 \KM)Y^\l)Y^ ) {2) 



for any tensors Y {kl) and Y (k2> 
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to calculate the matrix element (ip a ; F\Xj T \ij) b ; F). For the case when X is a one-body 
operator, the result trivially reduces to a single-particle matrix element 

(Va|Xgg|^ 6 ) = (j b jbz, JJ Z \jajaz)(t b t bz ,TT z \t a t az }(a\\X^\\b) (18) 

where a denotes quantum numbers of states tp a without third components, a = (n a , l a ,jaj s a, t a ). 
This expression serves to define our conventions for reduced matrix elements. For a 2-body 
operator, X, the reduced matrix element, (a; F\\X^ J ' T ^\\b; F), is 

(^fix^'Pj^f) = WaMxgpjihMta 

= (jbjbz, JJz\jajaz)(tbhz: TT z \t a t az ) (a; F\\X^ \\b; F) 

with 

(a;F\\X^\\b;F) = ]T E E D UcJb3h, Jlh, J; F)D(t a t b t h ,T 1 T 2 ,T; F) 

h&F T1T2 J1J2 

x ((ah)J 1 T 1 \\X^\\(bh)J 2 T 2 ) AS , (19) 



v (XT)n, hh \ t rr\ _( a ^ J l J ^TiT lz \Xj J ^ %h\J2J2z,T 2 T 2z ) AS 
{(ahjJ^WX^ , \\{bh)J 2 T 2 ) AS = — — ' (20) 

\J 2 J 2z , J J z \JlJlz) {±2±2z, -L J- z\J-ll-lz] 



and 

D(j a JbJh,JlJ2,J;F)= £ £ 



(21) 

A similar definition follows for D(t a t b th,TiT 2 ,T; F). The subscript 'AS' implies the matrix 
element is antisymmetric under the permutation of the two-particle states. The difference 
between summation over ip^ G F and that over h £ F should be understood, 

E =E E E • 

ip h eF h&F j hz eF t hz eF 
When the neighbouring even-mass nucleus forms a closed shell, D becomes simpler : 



D(j a jbjh, J1J2, J] F) dosed = ^2j a + 1 vS~TT U ^ ajh 3 J<2 ' ^ 

where the subscript 'closed' means that jhz runs from —jh to jh, and U is a recoupling 
coefficient defined in [17 and is related to the 6-j symbol. Now in our case, J = 0, T = 
1, t a = t b = th = i, s a = s b = 8^ = i, and we consider only closed shells in j-space^, in 
which case we have a much simpler expression, 

2t7 -+- 1 

D{j a jbjh, hh, J = 0; i? ) c losed = •Wa *i„ji 2 ■ * + 1 ■ ( 23 ) 



* 3 This is not the case for j4=96, where there is one unclosed shell. Here we simply neglect contribution 
coming from this unclosed shell. 
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But in isospin space, we consider not only charge-symmetric N = Z cases but also neutron- 
rich cases, and therefore we should retain the original definition, 



DrfTuT^F) = D(^,T 1 T 2 ,T = 1;F) 

Ex ^ (2*02; ^thz\TlTi z ) (^tbz, \th z \J2J2z){T2T2z-, ^F Z \T\T\ Z 
t hz £FT lz T 2z 



{ 2 ^bz : lTjs I ^ 



2 "02/ 



With eqs.(^, 23, 24) here and eqs. (53-56) of ref.[|l7]], it is straightforward to calculate 
the reduced matrix elements of the axial-charge operator. 



3 Numerical Results 

It is convenient to represent the results in terms of a ratio 

s {rl> a \F\M x \ik\F) = {a\F\\M x \\b\F) 

-body -body 1 1 "/ 

where x = 'tree', In, 2n, 'loop' or '2-body'. Thus 5 X is the contribution of the x part 
relative to the 1-body contribution (impulse approximation): <5t re e is the ratio of the one- 
pion-exchange contribution evaluated in the soft-pion limit with no loop correction, 5\ n 
shows the form-factor effect of one-pion-exchange coming from the loop correction to the 
vertices, and gives the two-pion-exchange contribution. The 5\ oop stands for the total 
loop contribution compared to the impulse approximation: 

^loop = <5l7T + ^27T, (26) 

while ^2-body; which is also denoted by <5 mec , represents the total 2-body contribution: 

<5mec = ^2-body = ^trec + #loop- (27) 

Furthermore, we introduce r]\ oop as the ratio of the total loop contribution to the tree 
contribution: 

f?loop = ^loop/ ^tree (28) 

In Table |, we give two sets of parameters used in this calculation; "PMR" denotes 



parameters adopted in [1C], and "OBEPR" denotes those obtained from the Bonn OBEPR 
potential [p3j| . There is a small difference of 1.1 % in the pion mass between the two. There 
is also a larger difference in the pion-decay constant /„-: In PMR, an experimental value is 
adopted, while in using OBEPR, f n is taken from the Goldberger-Treiman relation and the 
fitted value of the nNN coupling constant of g^NN = 13.68353. We observe this difference 
in fa has a significant effect on <5 mec but a smaller effect on r/i oop 

For each parameter set, we have used three choices of short-range correlation func- 
tions; the first two are a simple step function with different cut-offs, gi(r) = 6(r — di), 
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Table 1: Parameters used in the present work 





r Ivirv 


Uorv-T rt 




93.0000 MeV 


86.3982 MeV 




139.5673 MeV 


138.0300 MeV 




938.2796 MeV 


938.2796 MeV 


9A 


1.26 


1.26 


di 


0.70693 fm 


0.71480 fm 


d 2 


0.5 fm 


0.5 fm 


<lc 


3.93 fm- 1 


3.93 fm -1 



Table 2: Contributions to S mec for various choices of parameters for the ls^ — ► 0p]/ 2 
transition in A=16 : g±(r) = 9{r—d\) with d\ = 0.5 m^ 1 ~ 0.71 fm, cj2{r) = 9{r—d2) with 
c?2 = 0.5 fm and 53 (r) = 1— jo{q c r) with q c = 3.93 fm -1 





PMR 


OBEPR 


9{r) 


9l(r) 


h(r) 


93 W 


9i (r) 


92 (r) 


93 (r) 


^tree 


0.407 


0.459 


0.480 


0.472 


0.535 


0.560 


Slw 


-0.027 


-0.056 


-0.080 


-0.039 


-0.080 


-0.112 


fan 


0.068 


0.120 


0.160 


0.091 


0.162 


0.216 


^loop 


0.041 


0.063 


0.080 


0.052 


0.081 


0.104 


^mec 


0.449 


0.522 


0.560 


0.524 


0.616 


0.663 


?7loop 


10.2 % 


13.8 % 


16.6 % 


11.0 % 


15.2 % 


18.5 % 



di = \m~ ~ 0.71 fm, d 2 = 0.5 fm. The third correlation function is g%(r) = 1— jo(Qc r ), 
where jo(x) = sin(x)/x and q c is fixed as q c = 3.93 fm -1 [18]. 

In Table ||, we show the contributions to S mec for the lsi/2 — * QP1/2 transition in ^4=16 
with various choices of parameters and correlation functions. Successive rows in the table 
give the soft-pion tree contribution, the l7r-loop correction, the 2-7r-loop correction, the total 
loop correction (sum of rows 4 and 5), and finally the total meson-exchange contribution 
(sum of rows 3, 4 and 5). The last row denoted by ?7ioop expresses the total loop correction 
relative to the 'tree' contribution. This correction is generally small ranging from 10% for 
a cut-off of d ~ 0.7 fm to 15% for d = 0.5 fm. 

A short-range (SR) correlation function is necessary for calculations such as these, 
but there is little guidance as to its choice. Inevitably therefore some model dependence 
is injected here representing short-range phenomenology that is under poor control. For 
long-range operators, such as those arising from pion-range tree graphs, this is not a serious 
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problem as the matrix elements are not strongly influenced by the choice of correlation 
function. For the operators from the shorter ranged loop graphs, however, we can get a 
factor of two difference depending on the choice. This is evident from the results given in 
Table 2. Any conclusions, therefore, have to be tempered by this reality. 

A SR correlation function is required in this work for two reasons, and a common 
function is used here to cover both. First, there are correlations in the nuclear physics 
many-body problem. Because the two-body operators in eq. (||) are expressed in terms 
of the relative separation of two nucleons, r = \r% — f^|, it is convenient if the shell-model 
wavefunctions are expressed in terms of similar co-ordinates. This is quite practicable if the 
shell-model Hamiltonian is the harmonic oscillator. Then the coefficients of transformation 
from a single-particle basis to relative and center-of-mass basis are known. However, there 
is one drawback to the harmonic oscillator Hamiltonian: its eigenfunctions are not the 
eigenfunctions of the nucleon-nucleon interaction. This is particularly important in the 
relative coordinate, where the nucleon-nucleon interaction is known to have a strong short- 
range repulsion that makes the relative wavefunction go rapidly to zero as r — > 0, more 
rapidly than given by uncorrelated oscillator functions. Thus to incorporate this piece of 
many-body physics it is quite common to modify two-body operators by multiplying them 
by a SR correlation function, g(r), where g(r) is some function that tends to zero as r —* 
and tends to unity for large r. 

Second, in ChPT, loop integrations introduce high-momentum components. Al- 
though these loop integrals have been regularized and are finite, there are numerical difficul- 
ties in the Fourier transform to coordinate space. One way to deal with this is to impose a 
high-momentum cut-off. Another way is to impose a SR correlation function in co-ordinate 
space that would kill or ameliorate divergences at the origin. Consider a loop integral of 
order L. From the counting rules of ChPT and the nature of the loop integral, the most 
divergent piece has a form 



where /c is a non- negative integer depending on the process. In our case, k = 1 and L = 1. 
But, it is natural to assume that the correlation function is universal for all L. Thus, in 
order to ameliorate the divergence in all orders, we should have 



for any finite L. This is a rather severe restriction. In practice, because we are interested 
in one-loop accuracy, this restriction can be milder. In our case, we should have, at least, 



The step function 8(r — d) is the simplest SR correlation function and its underlying 
assumption for ChPT is that a nucleon has a hard core of radius d. It certainly satisfies all 





lim q(r) = or finite. 
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Table 3: Contributions to <5 mec with PMR parameters and with g(r) = gi(r) 





A=1Q 


A=40 


^4=48 


yl=96 


^4=132 


A=208 


^4=208 


■fc /l\ /T T 7"\ 

n(j(MeVJ 


13.921 


11.021 


10.489 


8.635 


7.874 


6.883 


6.883 


initial 


!Sl/2 


1P3/2 


!P3/2 


2si/ 2 


it 7/2 


igg/2 


2P1/2 


final 


fin 

U Pl/2 


Ua 3 / 2 


Oa-3/2 


-LPl/2 


U S7/2 


on 9 / 2 


On , 


°tree 


U.4U f 


U.4o0 


U.02o 


(\ AA1 
U.441 


U.004 


U.ooo 


U.4D2 




-0.027 


-0.031 


-0.032 


-0.028 


-0.034 


-0.035 


-0.026 


^2vr 


0.068 


0.077 


0.073 


0.069 


0.080 


0.084 


0.066 


^loop 


0.041 


0.046 


0.040 


0.041 


0.046 


0.049 


0.040 


<5mec 


0.449 


0.481 


0.564 


0.481 


0.600 


0.602 


0.502 


Vloop 


10.2 % 


10.5 % 


7.7 % 


9.3 % 


8.3 % 


8.8 % 


8.6 % 



the required properties stated above, it is easy to implement and it has a very clear physical 
interpretation. Similarly a step function has been used in the nuclear many-body problem as 
well [ 20 1 , where the interpretation of d now relates to the range of the short-range repulsion 



in the nucleon-nucleon interaction. 

There are other choices of SR correlation functions in use in nuclear physics. For 
example, Brown et al. [18] choose a form 1— jo(q c r), with q c = 3.93fm _1 , where q c has been 



adjusted to reproduce the dominant Fourier components of a realistic two-body correlation 
function calculated with the Reid soft-core potential. Note that the value of q c matches the 
inverse of the Compton wavelength for the w-meson. Thus this form of correlation function 



is particularly appropriate for any potential, such as the Bonn potential |23j, where the 
short-range repulsion is generated by 'enhanced' meson-nucleon coupling for the w-meson. 
The Bonn potential is 'softer' than the historic hard-core potentials and this correlation 
function is likewise softer than the form 6(r — d). If we define an effective cut-off, d e Q, for 
93(f) = 1 — jo(lcf), such that the correlation function 6(r — d e g) gives approximately the 
same result as 53 (r), then we obtain 

d eS ~ 0.4 fm. (29) 

This value is smaller than the choices of d\ — 0.7 fm and d<i = 0.5 fm used with gi(r) and 
92(f), respectively. 

The problem with the correlation function 1 — jo(q c f) is that it does not meet the 
requirements from ChPT of ameliorating the divergences at the origin. In particular it 
could not be used for two-loop or higher calculations. Thus in the following discussions we 
will give less weight to calculations using 53 (r). 

In Table [3|, we display the contributions to <5 mec for various odd-mass nuclei of closed- 
shell-plus-one configuration, where the mass number of the closed shell, A, is indicated. 
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The purpose is to investigate the mass dependence of the two-body contributions. Here, the 
PMR parameter set with gi(r) is adopted. In all cases, the l7r-loop and 2-7r-loop corrections 
are small and of opposite sign so that the resultant 5i 00 p is small. Further, when expressed 
relative to the 'tree' contribution, the 77i 00 p shows little mass or state dependence. We will 
discuss mass dependence further in the next section. 

In Table [|, we present our results for 5 mec and i]\ oop for various single-particle transi- 
tions in a number of different nuclei ranging from light to heavy, with the PMR parameters 
and the g\{r) correlation function. 
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Table 4: Calculations of S mec and r/i oop for a range of transitions in various nuclei with 
PMR parameters and including the short-range correlation function, g(r) = gi(r) = 6(r— 
with d\ = 0.5 m^ 1 ~ 0.7 fm 



A= 


=16, ftw = 13.921 


A=40, ?kj 


= 11.021 


Transition 


^mec 


??loop 


Transition 


^mec 


??loop 


0d 3 / 2 — ► 


0P3/2 


0.351 


9.2 % 


Of 5/2 — > 0d 5 / 2 


0.350 


9.4 % 


ls l/2 


OPl/2 


0.449 


10.2 % 


lp-J/9 — ► Odo/n 
i o/ z o/ z 


0.481 


10.5 % 


Pl/2 


Osi/2 


0.508 


9.3 % 


!Pl/2 ls l/2 


0.389 


10.1 % 










0d 3 / 2 -> 0p 3/2 


0.475 


9.4 % 










lSl/2 0p 1/2 


0.536 


10.2 % 


A= 


=48, hoj = 10.489 


A=96, hu. 


= 8.635 


w 0/ z 


0d 5/2 


0.409 


8.4 % 


°g7/2 0f 7/2 


0.386 


8.8 % 


IDo /o — > 


0d 3/2 


0.564 


7.7 % 


ld 5/2 -► 0f 5/2 


0.568 


8.4 % 


lPi /o —> 

11/ Z 


lSl/2 


0.434 


8.8 % 


ld 3 / 2 -> lp 3 / 2 


0.400 


9.3 % 


Odq/o — ► 
o/ Z 


°P 3 /2 


0.507 


9.4 % 


2Si/2 — lPl/2 


0.481 


9.3 % 


lSl /9 — ► 
1/ z 


OPl/2 


0.600 


8.4 % 


0f 5/2 -> 0d 5/2 


0.468 


9.4 % 










!p 3 /2 -»• 0d 3/2 


0.593 


8.7 % 










!Pl/2 ls l/2 


0.502 


9.2 % 




=132, 


fiu = 7.874 


A=208, ftu; = 6.883 


0h 9/2 -> 


0g 9 /2 


0.400 


7.8 % 


Oin/2 0h 11/2 


0.390 


8.2 % 


l f 7/2 ->■ 


0g 7 /2 


0.600 


8.3 % 


lgQ/2 ~> 0h 9/2 


0.602 


8.8 % 


lf 5 /2 - 


ld 5 /2 


0.406 


7.7% 


lg7/2 lf 7/2 


0.395 


8.1 % 


2 P3/2 ~ ► 


ld 3 /2 


0.517 


8.3 % 


2d 5 / 2 — »■ lf5/2 


0.530 


8.8 % 


2 Pl/2 


2si/ 2 


0.437 


8.1 % 


2d 3 / 2 — ► 2p 3//2 


0.418 


8.5 % 


0g 7 /2 -» 


0f 7 /2 


0.482 


8.3 % 


3Si/2 — > 2p x /2 


0.476 


8.8 % 


ld 5/2 -» 


0f 5 /2 


0.610 


8.9 % 


0h 9/2 -> 0g 9/2 


0.461 


8.5 % 


ld 3 /2 


1P3/2 


0.501 


8.1 % 


lf 7 / 2 — > 0g 7 / 2 


0.608 


9.3 % 


2s 1/2 -► 


!Pl/2 


0.559 


8.4 % 


l f 5/2 — ^ ld 5 / 2 


0.475 


8.4 % 










2p 3 / 2 !d 3 /2 


0.561 


8.9 % 










2Pi/2 - >■ 2si/ 2 


0.502 


8.6 % 
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From these tables, two conclusions emerge: 

The loop correction (??i 00 p) is around 10 % with some dependence on the choice of SR 
correlation function. This indicates quite a small correction, confirming the dominance 



of the soft-pion-tree graph and the conclusions made in [ 10 1 , where a simple-minded 



Fermi-gas model was adopted. The chiral filtering conjecture appears to hold. 

• The loop correction (^i 00 p) is essentially nuclear-mass and state independent. This 
is easily understood in that the tensor structure is common to the 'tree', In and 
part of the 2ir operators displayed in eq.(||). These operators only differ in their 
radial functions. Furthermore in light nuclei with LS closed shells, the other part 
of the 2-7T operator also has the same tensor structure because the matrix element of 
(<?i + <?2) {j\ x T2) is the same as that of (<ti x 02) (ti +T2) in this case. From Tables 3 
and 4, we observe that ryioop of neutron-rich nuclei is smaller than that of N = Z 
nuclei by 10 or 20 %, which may be viewed as the difference in the matrix element of 
the two operators for neutron-excess orbitals. 

4 Discussion 

We should like to conclude with a more detailed discussion of the mass dependence 
in the results. To this end it is convenient to define a ratio, r, where 

<5 mec ( 208 Pb) . . 

r =W (30) 

Here 5 mec is the correction expressed as a fraction of the 1-body impulse- approximation 
matrix element. For light nuclei in the vicinity of the closed shell, A=16, we consider the 
transition, lsi/2 — * Opi/2; that dominates the spectroscopy in this mass region, while for 
heavy nuclei in the vicinity of the closed shell, ^4=208, we consider the lgg/ 2 — > 0h g / 2 
transition. This ratio was introduced by Towner Q, because it is rather insensitive to the 
choice of SR correlation functions and the parameters used. Since this choice represents 
one of the biggest uncertainties in the present work, the ratio, r, is an attractive quantity 
to discuss. 

There is some experimental information on the ratio, r. It derives principally from ex- 
perimental data on first-forbidden beta decays, as analyzed in the shell model by Warburton 
et al. Q. The method is to compute beta-decay matrix elements in impulse approximation 
with the best available shell-model wavefunctions, and then allow the matrix element of the 
time-like part of the axial current to be multiplied by an enhancement factor, e mec - Note 
that 

e __ (-Ml-body) + (-Ms-body) _ j. , . 

tmec — - — 1 ~r O me c- v" 1 / 

(M l-body) 

A value for e mec is obtained from a fit between experiment and calculation over a number 
of transitions in the mass region under study. The principal difficulty in this analysis is the 
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Table 5: Table of r-values. The asterisks on the entries are a reminder of the special meaning 
of <5mec) explained in the text 





<W( 16 o) 


<W( 208 Pb) 


r 


comments 






tree 
loop 
mec 


0.407 
0.041 
0.449 


0.553 
0.049 
0.602 


1.36 
1.20 
1.34 


table 3 line 1, present work 
table 3 line 4, present work 
table 3 line 5, present work 


tree 
correction 
mec 


0.560 
0.031 
0.591 


0.759 
0.055 
0.814 


1.36 
1.77 
1.38 


table 2, Towner || 
table 2, Towner || 
table 2, Towner || 


c 
1] 
t 




BR scaling I 
BR scaling II 


0.567* 
0.595* 


0.885* 
0.885* 


1.56 
1.49 


Kubodera-Rho [0, 


25 


1 


Experiment 
Experiment 


0.61(3) 
0.61(3) 


0.79(4) 
0.91(5) 


1.30(9) 
1.49(10) 


weak tensor force, Warburton || 
strong tensor force, Warburton || 



inevitable truncation required in the model space used in the particular shell-model calcu- 
lation. Thus Warburton computes a correction for model-space truncations to first order in 
perturbation theory. This correction, however, is dependent on the choice of residual inter- 
actions used in the calculation. In particular, it depends quite sensitively on the strength 
of its tensor component. For weak tensor forces, such as obtained in the Bonn interaction 



[23], the value of r is 1.30 ± 0.09, while for strong tensor forces, such as obtained with the 
Paris potential [24], the value is 1.49 ± 0.10. In either case there is more enhancement of 
the axial-charge matrix element required in the lead region than in the oxygen region. 

In Table || we list some calculated r-values from the present work and other sources 
|, 0]. The first row gives just the one-pion tree-graph contribution, where a value of 
r = 1.36 indicates that most of the experimental mass dependence is accommodated by the 
shell model with just the soft-pion 2-body operator. The mass dependence arises because 
the number of core orbitals being summed over in eq.(|ll]) is changing from light to heavy 
nuclei and the oscillator frequency parameter is reducing to reflect the increase in nuclear 
size. One can easily get similar effects in Fermi gas models. For example, Delorme @ was 
the first to write down the expression for <5tree for a nucleon in a Fermi gas 



m N k F 

'tree 



1 + 2x - 2x(l + x) ml 1 + - 



(32) 



where x = m^/(4/cp) and k F is the Fermi momentum. In Fermi gas models, A; F is related to 
the nuclear density, p = 2/c|/(37r 2 ). If it is assumed that a valence nucleon in a light nucleus 
such as oxygen experiences only one half the nuclear matter density, then the appropriate 
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value of k F to use in eq.(|3^) is k F = 1.08 fm . For a valence nucleon in a heavy nucleus 
such as lead the appropriate nuclear matter value of k F = 1.36 fm -1 would be used. Then, 
with the PMR parameter set, <5 mec ( 208 Pb) = 0.54 and <5 mec ( 16 0) = 0.39 and an r-value of 
r = 1.4 is trivially obtained.^ 3 



A similar idea has been discussed by Kubodera and Rho |11J| . They argue from [21] 
that incorporation of approximate chiral and scale invariances of QCD leads to a chiral 
Lagrangian of low-energy hadrons in which the pion-decay constant and the hadron masses 
scale universally as a function of the matter density p according to: 

(33) 

m N m a m p f w 
while the axial-coupling constant does not scale, 

^ sa 1. (34) 
9A 

Here the asterisk refers to a value in a nuclear medium as opposed to the free hadron value. 
Note that the pion-decay constant, f n , and the heavy meson masses all scale according to 
the same function. However the mass of the pion, being a Goldstone boson, is assumed 
not to scale, m* ~ m n . These assumptions are referred to as Brown- Rho (BR) scaling. 
But there is a subtlety explained in detail in Rho's lecture note p5|]: the Gamow- Teller 
coupling constant §a scales in medium not due to the above BR scaling but due to the 
short-range interactions between baryons; as a result, the constant g\l f* associated with a 
pion exchange remains constant, 

% « f (35) 

J IT J IT 

These considerations lead to 

(A°(l-body)>* = Hp)- 1 (^°(l-body)), (36) 

and 

(^°(mec))* = Hp)' 1 (A°(mec)), (37) 
where the asterisked quantities should be evaluated with the scaled parameters defined in 



eqs.(|33|, p4 |35| ), while un-asterisked quantities should be evaluated with the parameters in 
free space, given in Table 1. Once the 1-body operator becomes p-dependent, we need to 
elaborate on the definition of e mec [eq.pl])] by specifying what 1-body matrix element is 
used in the denominator. To be consistent with the way in which the empirical e mec was 
deduced by Warburton et al. || |(|, we must adopt the definition 

_ (.4° (1-body))* + (^°(mec))* 



(A°(l-body)> 



(38) 



# 4 With the OBEPR parameter set, these values become <5 mcc ( 20S Pb) = 0.63, <5 mcc ( 16 0) = 0.45 and r = 1.4. 
The r-value is insensitive to the parameter set used, which can be understood by noting that the r-value 
depends only on the value of m w . 
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The S mcc which should appear in eq. 
by 5mec b is numerically defined by 



\) and which for the sake of clarity is denoted here 



rWarb 



Meanwhile, eqs. 



and 



imply that 



Hp)' 



(A (mec)) 
<A°(l-body)) 



(39) 



(40) 



Reinterpreting S mec of eq.pq) in the present context, we identify 

(4°(mec)> _ x 



<A°(l-body)> 



(41) 



Then 



srWarb 



(1 + S mec ) - 1. (42) 

We assume /)( 208 Pb) = po with po — 0.17 fm~ 3 being the normal nuclear matter density. 
There is latitude in choosing p( ie O), and we consider here two cases. In Case I we assume, 
as previously, that p( 16 0) = pq/2. To illustrate the sensitivity of the results to p, we 
also consider Case II, in which p( 16 0) = 0.6/3o- A typical choice of $ is <£(/?) = 1 — 
0.15 (p/po). Then, using (^mec given in Table we obtain, for Case I, ^rn^k(^C') — 0.567 
and 5^ c rb ( 208 Pb) = 0.885, which leads to r = 1.56. The corresponding numbers for Case II 
are: C^c b ( 16 0) = 0.595, Cec b ( 208p b) = 0.885, and r = 1.49. These results are listed in 
Table ||| under the heading 'BR scaling'. We note that in this approach the scaling factor 
&(p) is the principal mechanism for explaining the observed mass dependence. 

If we repeat the same procedure using <5 tree instead of <5 mec , the resulting r-value will 
be somewhat larger; r tre e = 1-59 for Case I, and r tTee = 1.51 for Case II. As can be seen 
in Table || the corrections coming from pion loops evaluated in ChPT are almost nuclear 
mass independent, and this explains why the sum of tree graphs plus corrections leads to a 
smaller r-value than the tree graphs alone. This feature is in contrast to the results obtained 
by Towner || in a quite different approach. There the tree graphs are not evaluated in the 
soft-pion limit but are evaluated with full momentum dependence retained, including vertex 
form factors^. This calculation, called the 'hard-pion approach', reduces the value of the 



tree graph significantly. The reduction, however, is largely compensated by heavy-meson 
pair graphs, which also give a sizeable contribution as pointed out by Kirchbach, Riska and 
Tsushima 0. Thus the correction in this approach is computed from 'hard pions' plus 
'heavy mesons' minus 'soft pions', and is given in line 5 of Table ||. This correction has 
significant mass dependence, such that the r-value of the sum of soft-pion tree graph plus 
corrections increases relative to soft pions alone. 



#5 We remark that form factors built into the Lagrangian have no place in ChPT; the momentum depen- 
dence of vertices in ChPT arises from radiative corrections, as one calculates higher-order diagrams in chiral 
expansion. 
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These two viewpoints, however, can be reconciled. When heavy-meson pair graphs are 
explicitly computed, the largest contributions come from a and u mesons. Mathematically 
the u-meson contribution takes the form of the 1-body impulse approximation with the 
nucleon mass replaced by an effective mass. This was shown by Delorme and Towner [ 13 1 in 
a different context. Recently Birse [26] finds a similar result in the non-topological soliton 
model for a nucleon embedded in mean scalar and vector fields. This implies that the 
phenomenological cr-meson and the BR scaling play the same role, although their chiral 
properties are quite different. At this place, it should be noted that there is no way to 
introduce a cr-field in a theory where chiral symmetry is non-linearly realized. The only 
possible way is to define a chiral-scalar (and of course Lorentz scalar) field as Brown and 
Rho did in their BR scaling [21|. Therefore to make comparisons with ref.pl, BR scaling 
has to be added to the ChPT results. 

In conclusion, we find that corrections to the soft-pion tree graph computed from 1-7T 
and 2-7r loop graphs in ChPT are small and around 10 %. This lends further support to 



the chiral filtering conjecture M, 22]. Second, the mass dependence or density dependence 



evident in the analysis of the experimental data by Warburton et al. Q has a variety of 
interlinking explanations ranging from the trivial mass dependence inherent in the shell 
model to the more fundamental role of heavy mesons or BR scaling. 
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